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Abstract 

We analyze constraints on the parameters of the R p violating su- 
persymmetry ($p3USY) which can be extracted from non-observation 
of the neutrinoless nuclear double beta decay (Qu{3{3) at a given half- 
life lower bound. Our analysis covers a large class of phenomenolog- 
ically viable ^ P SUSY models. We introduce special characteristics: 
the SUSY sensitivity of a /3/3-decaying isotope and the SUSY reach 
of a Oi//3/3-experiment. The former provides a physical criterion for a 
selection of the most promising isotopes for SUSY searches and the 
latter gives a measure of success for a Oi//3/3-experiment in exploring 
the ^pSUSY parameter space. On this basis we discuss prospects for 
exploration of supersymmetry in various Oz^/3/3-experiments. 
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Observation of the neutrinoless double beta (0^/3/3) decay %Y — > + 
2e~ of some nucleus \ Y would be an unambiguous signal of the physics be- 
yond the standard model (SM) of electro-weak interactions [|TJ, 0. This 
process is forbidden within the SM since it violates lepton number (L) con- 
servation. 

During a long time the 0^/3/?-decay has been attracting great theoreti- 
cal and experimental efforts. A number of experimental collaborations are 
involved in searching for this exotic phenomenon (see for instance 0). Un- 
fortunately, there is not yet any evidence for the 0z//?/3-decay. Nevertheless, 
impressive progress has been achieved in establishing a lower bound on the 
half-life of various isotopes and further advance in this direction is expected 
in the near future. 

The question is: what sort of new information on the physics beyond the 
SM we are provided with this experimental progress. 

It was a common practice to answer this question in terms of the Ma- 
jorana neutrino mass. This implies that the Oz//3/3-decay proceeds via the 
conventional mechanism based on Majorana neutrino exchange between the 
decaying neutrons, as presented in Fig. la. The measure of experimental suc- 
cess in this case is a reduction of the upper bound on the effective neutrino 
mass (m v ). The best result to date is (m u ) < 0.65 eV (90% C.L.) [|]. Other 
experiments are also nearly penetrating the sub-eV neutrino mass range. 

Quite recently it has been realized that the conventional neutrino ex- 
change mechanism is not the only possible one. Modern particle physics 
offers a mechanism based on the supersymmetric (SUSY) interactions 0-0. 
In a certain sense this mechanism is more interesting than the conventional 
one since it allows the /3/9-decay experiments in the exciting field of super- 
symmetry. SUSY interactions, inducing the Oz//5/3-decay, can contribute to 
many other processes. Therefore, one can compare information about these 
interactions obtained from different experiments including Oz//5/3-decay ex- 
periments. It has been shown in || that Oz//5/3-decay experiments are more 
sensitive to certain SUSY manifestations than the other running and forth- 
coming accelerator and non-accelerator experiments. 

In this note we give a quite general parametrization of the SUSY effect in 
the Oi//3/3-decay and formulate criterion for the success of a O^/5/3-experiment 
in terms of an upper bound on the effective SUSY parameter. The latter is 
a certain combination of the fundamental SUSY parameters. We introduce 
the sensitivity of an isotope to the SUSY manifestation and the SUSY reach 
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of a given O^/3/3-experiment characterizing the depth of its penetrating to 
the SUSY realm. The latter is represented by the unexplored part of the 
SUSY model parameter space. This provides us with a convenient basis 
for estimating prospects of the exploration of the SUSY in various 0u(3/3- 
experiments. 

We start with a short theoretical introduction. The 0i//9/9-decay is allowed 
within a special class of the SUSY models admitting B-L violation (B and 
L are baryon and lepton quantum numbers). These are presently popular 
SUSY models with R-parity violation (R p = ^—\^ B + L + 2S ( where S is the 
spin). The I]l, p terms can be introduced explicitly into the superpotential as 

= XijkLiLjEk + \ ijk LiQjD k + X' ijk UiDjD k . (1) 

The indices k stand for generations. L, Q denote lepton quark doublet 
superfields, while E, U, D correspond to lepton and up, down quark singlet 
superfields. The first two terms in (|l|) lead to lepton number violation, while 
the last one violates the baryon number. For the 0^/3/?-decay only the A and 
A' type couplings are of relevance. The presence of lepton number violat- 
ing interactions in the Jftp SUSY model allows one to construct the SUSY 
mechanism of the 0^/9/9-decay ||-|§. 

A complete analysis of this mechanism within the Ijtp minimal supersym- 
metric standard model (_^ P MSSM) was carried out in ||. On rather general 
grounds it was also shown || that the dominant contribution to this mech- 
anism within a phenomenologically viable ^? P SUSY model comes from the 
gluino g exchange diagram presented in Fig. lb. This contribution does not 
depend on specific details of a ^SUSY model such as the neutralino content, 
the mixing coefficients, etc. Thus, for any such model we can write down a 
dominant term of the effective quark-electron Lagrangian in the same form 
as for the $ P MSSM derived in 



£e//0) = i^zr ■ Vsusy 



G 2 r i 

JpJp + JsJs — ' Jthv 



2m p 

+ subdominant terms. 



(e(l + 75 )e c ) (2) 



The quark currents are defined as Jp = u ar y 5 d a , J s = u a d a , J^ v = u a cr^ lJ (l + 
j5)d a . The effective parameter t/susy accumulates certain fundamental pa- 
rameters of the ^SUSY model. In the particular case of the I/i p MSSM it 
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has the form 



\' 2 



m. 



9 G\m% m~a 

do 



(3) 



Here m g is the gluino mass and m - j l are masses of the superpartners of 
ul and quarks; a s is the strong coupling constant. Starting from the 
Lagrangian (0) describing the basic quark-lepton transition d + d — > u + 
u + 2e~, which triggers the 0^/3/3-decay, one can derive the following half-life 
formula J7fl 
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The factor in brackets implies the normalization of M.q as in ||. = 
0.85 GeV is the momentum scale of the nucleon dipole form factor F(q 2 ) = 
(1 + q 2 /m A )~ 2 used in calculation of the nuclear matrix element Aig. The 
leptonic phase space integral G i is defined as 



G, 



01 



(G F cose c fA/fv)*m* e 
64vr 5 ^(m eJ R ) 2 ln2 



dVtid^bm 



(5) 



and tabulated in for various nuclei. Values of G i for nuclei analyzed in 
the present letter are given in the Table. Here, Rq is the nuclear radius and 
Ja = 1.261, fy = 1. The kinematical factor b m accounts for the Coulomb 
distortion of the electron waves ||. 

The nuclear matrix element Aig for the SUSY contribution to the Qv/3j3- 
decay can be written in the form independent of the nuclear wave function 



Mr, ^<0f || £t-W } ( — J \ocPf n (x a ) + a$F A {z A )+ 



(6) 



+ °i ■ <Jj (uaF n (x a ) + a { A ] F A (x A )) + a r {3(^ • r ij ){a j ■ r^)- 
- a i -a J }F 5 (x A )]||0+). 

The following notations are used 

r ij = ( r i ~ r j ) > r «i 



| Fiji, TijHijy "^A^iji 
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where ~r*i is the coordinate of the ith nucleoli. The nucleon structure pa- 
rameters ay 1 = 0.15, ay 1 = 7.5, a^' = —1.2, a^' = 0.28, ax = 1.3 were 
calculated in [7J. 

Three different structure functions Fi are given by the integrals over the 
momentum q transferred between two decaying nucleons and have the form 

Fn(x) = ^(3+3x+x 2 ), F 4 (x) = ^!(3+3x-x 2 ), F 5 (x) = ^1. (7) 

Expression (|6|) allows numerical calculation of the matrix element A4q within 
any nuclear structure model. So far such calculations have been carried out 
only within the pn-QRPA (proton-neutron Quasiparticle Random Phase Ap- 
proximation) in ||. It is well- motivated and reliable approach ]10 success- 
fully describing the 2z//3/?-decay recently observed in several experiments. 
Numerical values of M.q for several experimentally most interesting isotopes 
are given in the Table. We point out that this matrix element drastically 
differs in structure from the well known Qvf3f3 nuclear matrix elements of the 
Majorana neutrino exchange mechanism. 

Non-observation of the 0z//3/3-decay at the half-life level Ty£ leads to the 
following constraint on the ^ P SUSY parameter rj SUSY 



10" 7 

J 1/2 > I l/2 * VSUSY S 



C(Y)\ 



1 -^S P ^ = lO- 7 -e(Y,ex P r\ (8) 

J l/2 



where Tfj 2 is the theoretical expression given in (|j). We have introduced the 
following characteristics of a /3/3-decaying isotope Y and a particular Ois/3/3- 
experiment giving a quantitative basis for assessing its ability to search for 
a SUSY signal. These are the SUSY sensitivity of an isotope Y 



OX) = 10 5 M^G 01 (9) 
and the SUSY reach of the 0^/?/3-experiment with the isotope Y 



\ 10 /4 years 



rpexp 

1 2 (io) 



The former is an intrinsic characteristic of an isotope Y depending only on 
the matrix element and the phase space factor. The latter characterizes the 
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experimental set-up as a whole. TfJ% is the lower half-life bound reached 
with this set-up. This is a time dependent characteristic for improving ex- 
periments. 

The large numerical values of the SUSY sensitivity ( in (|9|) correspond 
to those isotopes within the group of /3/3-decaying nuclei which are most 
promising candidates for searching for SUSY in O^/3/3-decay experiments. 
Using general formulas in ©-(0) one can calculate ( in the framework of any 
nuclear structure model. Its numerical values calculated in the pn-QRPA 
are presented in the Table. As is seen, the most sensitive isotope is 150 Nd, 
then follows 100 Mo. Unfortunately, we are not able to make any conclusion 
concerning 48 Ca, which is also an experimentally interesting isotope. This is 
because the pn-QRPA does not provide reliable results for this light nucleus. 
It would be important to reconstruct a whole picture showing potential abil- 
ities of all /3/3-decaying nuclei from the point of view of their sensitivity to 
the SUSY signal. Therefore, calculation of the SUSY nuclear matrix element 
M.q for 48 Ca within a certain nuclear model is demanded. It could be done, 
for instance, on the basis of the nuclear shell. This approach allows exact 
nuclear structure calculations for 48 Ca fi~lfl , [ |i~2]1 . 

Of course, the SUSY sensitivity ( cannot be the only criterion for selecting 
an isotope for the O^/3/3-experiment. Other microscopic and macroscopic 
properties of the isotope are also important for building a O^/3/3-detector. 

The SUSY reach parameter e(Y, exp) of the O^/3/3-experiment takes into 
account both the SUSY sensitivity of the isotope Y and specific experimen- 
tal conditions determining the reach in the half-life limit. The current and 
the near future experimental situation in terms of the accessible half-life and 
the SUSY reach parameter e is presented in Fig. 2a,b. For completeness, 
in Fig. 2c we also present a situation with the experimental reach in the 
inverse neutrino mass. Fig. 2c is taken from [[|. Figs. 2 show the run- 
ning /3/3-experiments: the Heidelberg- Moscow (H-M) 76 Ge experiment [|], 
the NEMO 100 Mo experiment [O, and other 0^/9/9-experiments that also 
crossed the half-life barrier of 10 21 years with the isotopes 48 Ca fli~4 |; 82 Se [|15 
100 Mo 



TBI 



Cd |7|; 130 Te 0; 136 Xe (|TJ; 150 Nd [f?tj. We also included in 



Fig. 2 the forthcoming or partly operating experiments: the Gottard 136 Xe 
TPC experiment [21|, the 130 Te cryogenic experiment 



T|, the ELEGANT 



experiment with 64 g of Ca |22], the experiment with 1.6 kg of Nd using 



an improved UCI TPC [20 



In Fig. 2a two broken lines with black squares and triangles correspond 
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to the two fixed values of the SUSY reach e( 76 Ge, H-M)P resent = 5.3 (lower 

line) and e( 100 Mo, NEMO) future = 22.2 (upper line). Experiments reaching 
these lines provide the same constraints on t/susy in © as the Heidelberg- 
Moscow experiment at present (lower line) and as expected from the NEMO 
experiment in the near future (upper line) respectively. 

These lines are interrupted on the left from the 76 Ge since we have no 
value of the SUSY nuclear matrix element M.q for 48 Ca. 

It is seen from Fig. 2 that the Heidelberg-Moscow 76 Ge detector is the 
best m v explorer at present and up to the year 2000. This experiment is also 
the best SUSY explorer at present. However, in the near future the NEMO 
experiment is able to take over the leadership in exploration of SUSY as 
follows from Fig. 2a. 

In conclusion, we stress that from the physical point of view the lower 
half-life bound itself is not a satisfactory characteristic of the 0i//3/3-decay ex- 
periment. An objective characteristic should be related to the fundamental 
parameters of the underlying physics. In this note we give the parameteri- 
zation of the SUSY contribution to the 0^/3/3-decay in terms of one effective 
parameter rjsusY- The parameterization is valid for a wide class of ^? P SUSY 
models and can be used for the quantitative presentation of the SUSY effect 
in the O^/3/3-decay experimental data. Together with the effective Majorana 
neutrino mass (m u ) the parameter tjsusy is an important physical charac- 
teristic of the 0^/9/9-decay experiment. 

Our results for the SUSY sensitivities ( of various /3/3-decaying isotopes 
and for the SUSY reach e of /?/9-decay detectors based on the tjsusy param- 
eterization can be used for planning the future experiments. 
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Figure Captions 



Fig.l (a) The conventional massive Majorana neutrino exchange mechanism 
and (b) the supersymmetric mechanism (the dominant contribution) of 
the neutrinoless double beta decay. 

Fig. 2 1995 experimental situation (gray bars) and expectations up to 2000 
(open bars), (a) half-life Tu 2 ; two broken lines correspond to the two 
fixed values of the SUSY reach e( 76 Ge, H-M)P resent = 5.3 (lower line) 
and e( 100 Mo, NEMO)^ ure = 22.2 (upper line); experiments reaching 
these lines provide the same constraints as the Heidelberg-Moscow ex- 
periment now (lower line) or as expected from the NEMO experiment 
in the near future (upper line), (b) the SUSY reach parameter e (see 
(|10D ) and (c) the inverse neutrino mass reach l/m u (taken from J|]). 
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Table: The SUSY sensitivity (, phase space factor G m and nuclear matrix 
element for the SUSY mechanism of the Oi//3/3-decay calculated within the 
pn-QRPA for several experimentally interesting isotopes. ( defines the con- 

( 24 V /2 

straint on the SUSY parameter r/susY < -rTrf Jo%/3 



Isotope 

Ay 


Nuclear Matrix 
Element 

M~ q 


Phase Space 
Factor 
Goi • 10 14 


SUSY 
Sensitivity 

c 


150 Nd 


416 


20.808 


19 


100 Mo 


328 


4.56316 


7 


130 Te 


262 


4.41596 


5.5 


82 Se 


253 


2.80583 


4.24 


116 Cd 


190 


4.92614 


4.2 


136 Xe 


143 


4.71196 


3.11 


76 Ge 


283 


0.635941 


2.26 


128 Te 


298 


0.181888 


1.27 
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